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Abstract 

We introduce a model for a growing random graph based on simultaneous reproduction 
of the vertices. The model can be thought of as a generalisation of the reproducing graphs 
of Southwell and Cannings and Bonato et al to allow for a random element, and there 
are three parameters, a, (3 and 7, which are the probabilities of edges appearing between 
different types of vertices. We show that as the probabilities associated with the model vary 
there are a number of phase transitions, in particular concerning the degree sequence. If 
(1 + a) {I + 7) < 1 then the degree distribution converges to a stationary distribution, which 
in most cases has an approximately power law tail with an index which depends on a and 7. 
If (1 + a)(l + 7) > 1 then the degree of a typical vertex grows to infinity, and the proportion 
of vertices having any fixed degree d tends to zero. We also give some results on the number 
of edges and on the spectral gap. 

AMS 2000 Subject Classification: Primary 05C82; secondary 60G99, 60J10 

Key words and phrases: reproducing graphs, random graphs, degree distribution, phase 

transition 



1 Introduction 



In this paper we introduce a new model for a growing random graph based on simultaneous 
reproduction of the vertices in the graph, with edges being formed between the new vertices and 
each other and between the new vertices and the existing ones according to a random mechanism 
conditioned on the pattern of edges between the vertices in the previously existing graph. The 
model is a generalisation of the models introduced by Southwell and Cannings [11, 12, 10] and 
the Iterated Local Transitivity (ILT) model of [2], introducing stochasticity, which causes the 



1 



regular structure found in the graphs of [11, 12, 10] to be lost, and which may make them more 
suitable for modelling in areas such as social networks; the authors of [2] particularly suggest 
their model as a model for online social networks, mentioning Facebook and Twitter among 
other examples. 

We will show that our model, which depends on three parameters a, f3 and 7, which are to 
be thought of as probabilities, exhibits a number of phase transitions as the parameters vary; 
for example for some values of the parameters we will show that the degree distribution of the 
graph converges to a limiting probability distribution, while for other choices of the parameters 
the degree of a randomly chosen (in an appropriate sense) vertex in G n can be shown to tend 
to infinity as n — > 00. We will also show that for certain choices of the parameter values the 
model exhibits a power-law-like decay of the degree distribution, which is a property reported 
for many "real world" networks, and is also associated with other random graph models such as 
preferential attachment. 

We start with a graph Go, and form a new graph G n+ \ by adding a "child" vertex for every 
vertex of G n . As in [11, 12, 10] we denote the vertices by binary strings, writing vO for the "child" 
of vertex v 6 V(G n ) and vl for the continuation of vertex v as a vertex of G n+ \. The edges of 
G n +i are then obtained according to the following mechanism. For each n define independent 
(of each other and of the random variables at other stages of the construction) Bernoulli random 
variables a$ v \ ~ Ber(a) for each unordered pair {u, v} of vertices of G n , ~ Ber((3) for 

each vertex in G n , c|"^ ~ Ber{^) for each ordered pair (u,v) of vertices of G n , and connect 
vertices as follows: 

(a) ul is connected to vl in G n+ \ if and only if u and v are connected in G n , that is existing 

edges are retained, and no further edges are formed between existing vertices. 

(b) uO is connected to ul in G n+ \ if and only if = 1, so each child is connected to its parent 

with probability f3. 

(c) uO is connected to vl in G n+ \ if and only if c|™^ = 1 and u and v are connected in G n , so 

each child is connected to each of its parent's neighbours with probability 7. 

(n) 

(d) uO is connected to vO in G n+ \ if and only of a^J v y = 1 and u and v are connected in G n , 

so each child is connected to each of its parent's neighbours' children with probability a. 

The models introduced in [11, 12, 10] have a, (3, 7 £ {0,1}, so are deterministic. Additionally 
the case where a = 0, /3 = 1, 7 = 1 is the ILT model, introduced in [2] as a model for online social 
networks. The ILT(p) model introduced in [2] as a stochastic generalisation of the ILT model 
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adds extra random edges between the child vertices without regard to whether their parents 
were connected, and thus cannot be seen as a special case of our model. In addition as defined 
in [2] the ILT(p) model always has at least the edges found in the basic ILT model, so is not 
suited to producing relatively sparse graphs. 

The model differs from duplication graphs, for example those considered in [3, 5], in that in 
those models only one vertex, chosen at random, duplicates at any one time step, whereas in 
the models considered here and in [2, 11, 12, 10] all vertices simultaneously duplicate. 

Our main results concern the degree distribution. We will deal with the cases where (3 = and 
(3 > separately, as the behaviour of the model when f3 = is potentially quite different, with 
large numbers of isolated vertices. 

Theorem 1. Let (3 = 0. Then, if (1 + 7) (a + 7) < 1, the probability that a randomly chosen 
vertex in the graph G n is isolated tends to 1 as n — > 00, and the proportion of vertices in the 
graph with degree zero tends to 1, almost surely. If (1 + 7) (a + 7) > 1, then the probability that 
a randomly chosen vertex in the graph G n is isolated converges to some value strictly less than 
1. 

Theorem 2. Assume (3 > 0, and letp^ be the proportion of vertices in G n with degree d. Then 



(a) // (1 + 7) (a + 7) < 1 there exists a random variable X such that p£ — > P(X = d) as 

n — > 00, almost surely. 

(b) Under the conditions of (a), the random variable X has a finite pth moment if (1 + 7)^ + 

(a + 7) p < 2, and does not have a finite pth moment if (1 + j) p + (a + 7) p > 2. 

(c) // (1 + 7) (a + 7) > 1 then p^ — >■ as n — >■ 00 , almost surely. 



Note that Theorem 2(b) implies that if (1 + 7) p + (a + 7) p = 2 the tail of the degree distribution 
is asymptotically close to a power law degree distribution with index given by — (p + 1), in the 
sense that qth moments exist for q < p but not for q > p. 

In [2], it is shown that the ILT model exhibits a "densification power law", which is defined 
to mean that, if E n is the number of edges of G n and V n the number of vertices, then E n is 
proportional to (V n ) a for some a G (1,2). The following result shows that our model exhibits a 
phase transition in this respect, with the transition occurring where 27 + a = 1. Note that in 
our model, as in the ILT model, V n = 2 n Vo for all n. 

Theorem 3. (a) If 27 + a > 1 then W n = ( 1+2 ^ a )n converges to a positive limit, so that the 
model has a densification power law as defined by [2] with exponent '"^^J" 1 ""- 1 . 
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(b) // 2 7 + a < 1 then 

En ^ V P 

2™ ^ 1 - 27 - a ' 

almost surely, as n — > 00, so £/ia£ i/ie number of edges grows at the same rate as the number 
of vertices 

(c) // 27 + a = 1 i/ten 

2 n n 2 ' 

almost surely, as n — >■ 00. 



Note that the combination of Theorems 3 and 2 implies that when 27+a > 1 but (l+7)(a+7) < 
1 the process exhibits both a densification power law in the sense of [2] and an approximately 
power law limit for the degree distribution. 

A further result in [2] on the ILT model concerns the spectral gap. They show that the normalised 
graph Laplacian C, as defined by Chung [4], of the ILT model has a large spectral radius, defined 
as max{|Ai — 1|, |A n _i — 1|}, where Ai is the second smallest eigenvalue (the smallest being Ao = 
for any graph) and A ra _i is the largest eigenvalue and thus that the graph has relatively poor 
expansion properties. The following results show that the same is also true for our model. We 
concentrate on the case /3 = 1, where the graphs are connected; otherwise Ai will be zero. The 
proofs use the Cheeger constant and its relationship to Ai, as defined in Chapter 2 of Chung [4]. 

Theorem 4. Let (3 = 1 and assume that Go is connected, so that G n will also be connected for 
alln. Let \i{G n ) be the smallest non-negative eigenvalue of the Laplacian of G n . Then 

(a) 1/27 + a < 1 then \\{G n ) — > as n —> 00. 

(b) //27 + a > 1 then there exists a (random) A, with A < 1 almost surely, such that 

limsup^^ Ai = A. 



2 Pictures 



This section shows a few examples of graphs of this type, which were generated using a script 
written with the igraph package, [6], in R. The first two pictures show examples with n = 7, 
a = 0, (3 = 1 and 7 = 0.2 and 0.49 respectively. 
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The next two pictures have a = 0, 7 = 0.366(? 
the second /3 = 0.8. 



and again n = 7. The first has j3 = 0.5, 
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3 Proofs of Theorems 



We start with a lemma on the conditional expectation and variance of the number of edges in 
G n , E n . This lemma will be useful for obtaining the mean of the stationary distribution of a 
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Markov chain which we will use to prove Theorems 1 and 2. We define T n to be the a-algebra 
generated by the graphs G m for m < n. 

Lemma 5. The conditional expectation and variance of E n+ \ satisfy 

¥.(E n+1 \F n ) = (1 + 2 7 + a)E n + 2 n /3V 
Var(£ n+1 |.F n ) = J E n (2 7 (l- 7 ) + a(l-a)) + 2 n y /3(l-/3) 

Proof. This follows from the fact that the E n+ \ can be written E n + E n+ \^ + E n+ \^ + E n+ \^ 
where E n+ i t i, E n+ i^ and -E n +i,3 are independent, E n+ \^ represents the edges between parents 
and children of their neighbours and, conditional on T n has a Bi(2E n ,~f) distribution, E n+ i^ 
represents the edges between children of neighbouring vertices and, conditional on T n has a 
Bi(E n , a) distribution, and £^+1,3 represents the edges between parents and their children and, 
conditional on T n has a Bi(V n , (3) distribution. As V n = 2 h Vq the result follows. □ 



3.1 Proof of Theorem 3 

We start with (a). The following approach is based on that in [1] for multitype branching 
processes, the idea being that the vertices and edges in the graph G n can be thought of as the 
two types in a population undergoing branching. However the resulting multitype branching 
process is not irreducible, so the results in [1] cannot be used directly. 

Given an edge in G m between vertices u and v, there will be an edge in G m+ i between ul and 
vi, and in addition there will be edges between ul and vO and vl and uO each with probability 
7 and an edge between uO and vO with probability a. We can consider these edges as offspring 
of the edge between u and v, and thus consider the set of edges in G n (for n > m) which are 
descendants of the edge between u and generation in a Galton- Watson branching 

process with offspring mean 1 + 2 7 + a, and where the extinction probability is zero and the 
number of offspring bounded. Treating the descendants of a given edge in G m as a subset of the 
edge set of G n , this shows that liminf ^ 1+ ^_ a ^ n is a positive random variable. 

Now define 

Vn i fl &n 

W = E 

(l + 2 7 + a) n 

ThenE(W n+ i|j;), so (W n ) n gpj is a non-negative martingale, and thus almost surely has a non- 
negative limit W. The above conclusion on liminf ^ l+ ^_ a ^„ shows that P(W = 0) = 0, giving 
the result. 



6 



For (b), Lemma 5 shows that E(E n ) = 1 _y_ ct 2" + o(2 n ) and 

Var(£ n ) = 2- 1 — ^ (2 7 (l- 7 ) + o (2")+a(l-a)) + 2«F /3(l-/3) + (l + 2 7 +a) 2 Var^x), 

1 — 2 7 — a 



which shows that 



max 



2™ J \\ \2' 4 
which implies the result via Chebyshev's inequality and the Borel-Cantelli Lemmas. 

For (c), an iterative use of Lemma 5 shows that E(E n ) = 2 n (^Eq + ^-n^ and Var(E' n+ i) 
2 n - 1 n/?y (27(l - 7) + a(l - a) + 4) + 0(2 n ). Then 



2 " n y \^n 2 2" / ' 
allowing the Chebyshev/Borel-Cantelli argument again. 

3.2 Proof of Theorem 4 



We consider some small m, and find the Cheeger constant of G m . By the definition in [4], this 
will be e(S m , S m ) / vol(S m ) for some S m C V(G m ), where for two subsets of the vertex set S 
and 5' e(S,S') is the number of edges between a vertex in 5 and one in S', and vol(5) is the 
sum of the degrees of vertices in S. (Note that vol(S') = 2e(S,S) + e(S,S.) Now consider the 
descendants of S rn in G n as a subset S n C F(G n ). Then the same arguments as in the proof 
of Theorem 3, applied to the subgraphs descending from S m and S m , show that if 27 + a < 1 
then the e(S n ,S n ) and e(S n ,S n ) both grow at rate 2™, in the sense that e ^f n) and e{S ^ Sn) 
converge almost surely to positive constants as n — > 00, and similarly if 27 + a = 1 e(S n ,S n ) 
and e(S n ,S n ) both grow at rate 2 n n, and if 27 + a > 1 e(S n ,S n ) and e(S n ,S n ) both grow at 
rate (1 + 27 + a) n . 

Next, again as in the proof of Theorem 3, (e(S n , S n )) n £^ forms a Galton- Watson branching 
process with mean of the offspring distribution 1 + 27 + a, and extinction probability zero, so 
e(S n , S n ) will grow at rate (1 + 27 + a) n . So for 27 + a > 1 min{v ol(s^li(S„)) ' which W the 
definition in [4] is greater than the Cheeger constant of G n , converges to a constant (less than 
1, as vo\(S n ) = 2e(S n ,S n ) + e(S n , S n )) and this constant bounds the lim sup of the Cheeger 
constant of G n above. 
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as n -> oo, and hence 
= O (^) , and again so is 

Hence by the Cheeger inequality (Lemma 2.1 and Theorem 2.2 of [4]), limsup,^^ \\{G n ) < 1 
almost surely in the case 27 + a > 1, and Xi(G n ) tends to zero in the case 27 + a < 1. 

3.3 Proofs of Theorems 1 and 2 

The proofs of Theorems 1 and 2 will rely on defining a certain Markov chain whose value X n 
represents the degree of a random vertex in the graph G n . We construct this by letting vq be 
a vertex of Go chosen uniformly at random, and then, using the binary string notation for the 
vertices described above, for n > 1 let v n = v n -±l with probability 1/2 and letting v n = v n -±0 
with probability 1/2. We then let X n be the degree of v n in G n . 

Then 

X n +1 = £,n+lX n + (1 — £n+l)W"n+l + Y n+ \ + Z n+ i, (1) 

where, conditional on G n , Y n+ \ ~ Bin(X n , r y), W n+ \ ~ Bin(X n , a), Z n+ \ ~ Bin(l,f3) and 
£n+i ~ Bin(l, 5), with all these variables being conditionally independent given G n . 

Here, ^ n+ i = 1 if our vertex in G ra+ i is a parent and if it is a child, W n+ ± represents child-child 
connections (so does not appear if £ n+ i = 1), Y n+ \ represents connections between a child and 
its parents' neighbours, and Z n+ i represents the connection between the child and its parent. 

As defined above, (X n ) ne m is a discrete time Markov chain on the natural numbers (including 
zero if /3 < 1). It is irreducible and aperiodic if /3 > 0, a < 1 and 7 < 1. (If (3 = then zero 
is an absorbing state, and if either a or 7 is 1 then X n is increasing in n and so the chain is 
certainly not irreducible, but otherwise P(X n+ i = l\F n ) is always positive.) 

Proposition 6. If 27 + a < 1 the distribution of X n converges in the Wasserstein-1 metric to 
a unique fixed point with finite mean . 

Proof. Note that if we have another random variable X n with a different distribution on No, we 
can apply (1) to it by defining, conditional on X n , Y n+ i ~ Bin(X n ,j) and W n+ \ ~ Bin(X n ,a) 
using the same set of Bernoulli trials as for Y n+ \ and W n +\ respectively, and letting 

X n +1 = £,n+lX n + (1 — Cn+l)Hn+l + Yn+1 + Zn+1, 



In the case where 2 7 + a < 1 . e &'? n) UQ ^ = O (( 1+2 J +a ) ) ->• 

' max(vol(5 n ),vol(5 n )) \\ 2 J J 

so is the Cheeger constant. Similarly if 27 + a = 1 , e |fg'f"\/ g : 

the Cheeger constant. 
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in which case E(|X n+ i — X n+ i||.F n ) = ^(27 + a + \)\X n — X n \, so that we have a contraction in 
the Wasserstein metric if 27 + a < 1. Hence in this case there is convergence in the Wasserstein-1 
metric of the degree distributions to a unique fixed point with finite mean. 

We can calculate the mean of this distribution by using Lemma 5: letting m = 2 we get 

E(4 n+1) ) = (2 7 + a + l)4 n) + 2 n /3v , 

where vq is the number of vertices in the initial graph, and solving this we find that the expected 
number of edg es in G n is 

/3v (2 n -(l + 2 7 + a)") 
1 - 2 7 - a 

so (as the number of vertices in G n is 2 n t>o) the expected average degree is 

/3(2"-(l + 2 7 + q)") 
2«-i(l-27-a) ' 

which converges to 1 _^_ a asru-oo. □ 



To go further than this we use Foster-Lyapunov techniques, as described in Meyn and Tweedie [8] 
in the more general case of an uncountable state space. The following lemma on the conditional 
moments of X n+ \ (including negative and fractional moments) will be useful. 

Lemma 7. Let p € R. Then as x — > oo ; 

E ((I±^±±^)V=*)-(l + 7)>, 

and 

E (( 1 + M -"+'^ + ^+' )V _ x) -» (« + 7F . 



Proof. In the case where p < this is a special case of Theorem 2.1 of Garcia and Palacios in 
[7]. When p > we can adapt the argument of Theorem 2.1 in [7] with the additional condition 
that the random variables are bounded above by a constant multiple of their mean, which is 
satisfied in this case. We assume that (A n ) n ^ is a sequence of positive random variables with 
associated sequences of constants jJL n and a n such that A n < kfi n with probability 1 for some k 
and for all n, and that the conditions of Theorem 2.1 in [7] hold. 



9 



For the upper bound, 

E<y*) = E(YVl {Yn< ^J+E(Yp {Yn ^ n+ ^ } ) 
with order of magnitude /4v 

For the lower bound (which is obvious by Jensen's inequality if p > 1), 

E(Y P ) > - tlifnXn >iln~ £), 

again with order of magnitude [)P n . □ 

Proposition 8. // (1 + 7) (a + 7) < 1, the Markov chain is positive recurrent, and thus the 
distribution of X n converges to a stationary distribution. 

Proof. This uses Theorem 11.0.1 of [8]. 

We choose p <G (0, 1) such that (1 + j) p + (a + j) p < 2. Because gj;((l + 7) p + (a + 7) p ) is negative 
at p = if log(l + 7) + log(a + 7) < 0, it will be possible to find such a p if (1 + j)(a + 7) < 1. 

We now let V(x) = x p . In [8], the drift AV(x) is defined as 

AV(x) = E(V(X n+1 ) - V{X n )\X n = x), 

and by Theorem 11.0.1 of [8] the chain will be positive recurrent if (for some V) AV(x) < —1 
for x large enough. Now 

E « +1 ,*„^> - f ( E (( 1 + I^ + ^)^) + «((»W^ + ^)V) 

< y((l+7) p + (a + 7)") + o(x") 
(by Lemma 7), 

so 

Ay(,)<^( (1 + ^ + (a + ^-l) +0 (^ 
which will be less than — 1 for x large enough, giving the result. 

□ 
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We now investigate the tail behaviour of the stationary distribution, in the case where Proposi- 
tion 8 shows one exists. 

Proposition 9. Let p > 0. If (1 + j) p + (a + j) p < 2, then a random variable X with the 
stationary distribution of the chain has finite pth moment K(X P ), and we have convergence of 
pth moments, E(X£) — >• E(X P ) as n — > oo. 



Proof. Again this uses a Foster-Lyapunov type technique, in this case Theorem 14.0.1 of [8] 
which states that if, for a given function / > 1, we can find V such that AV(x) < —f(x) for 
x large enough then / has a finite integral with respect to the stationary distribution and that 
E(f(X n )) converges to this integral. We will set f(x) = x p + 1. 



Let V(x) = kx p , where k is chosen so that 



•(i + ,r+ (a+T r_ 



Then, by Lemma 7, 



Ar(.r) < A-.,-'' ( ( 1+ T) P + (« + 7) P _ 1 j + o{ rl>y 



and so AV(x) < —f(x) for x large enough, giving the result. □ 

Proposition 10. Let p > 0. If (I + "f) p + (a + j) p > 2, then a random variable X with the 
stationary distribution of the chain does not have finite pth moment K(X P ). 



Proof. As Z n+ \ > 0, we have 

E(X p +l \X n = x) > °L (e (Ji + \ Xn = x ^ + E (JYL^L + ^±lj \ Xn = x)) , 

so by Lemma 7 

E(X£ +1 |X n = x) > ^ ((1 + 7 y + (a + 7 ) P ) + o(^). 

Hence the pth moment of X n tends to infinity as n — > oo, so by Theorem 14.0.1 of [8], again 
applied to f(x) = x p + 1, the stationary distribution cannot have a finite pth moment. □ 

Proposition 11. // (3 > and (1 + 7) (a + 7) > 1 the Markov chain is transient. 
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Proof. By Lemma 7, 



E((l + X n+i y\X n = x) (l + 7 ) p + (« + 7) p (t) s 
(1 + x)p 2 ' [) 

so we apply Theorem 8.0.2 (i) of [8] with V(x) = 1 — (1 + x) p for some p < such that 
(1 + 7) p + (a + 7) p < 2. With this choice of V, (2) shows that AV(x) > for x large enough, 
and as V is bounded and positive on the natural numbers Theorem 8.0.2 (i) of [8] gives the 
result. 

□ 



The case where (3 = is something of a special case as the chain is not irreducible. However we 
can show that when (l + 7)(a + 7) < 1 the probability that a randomly chosen vertex is isolated 
tends to 1, while there is positive probability that a randomly chosen vertex is not isolated when 
(l + 7)(a + 7) > 1. 

Proposition 12. // j3 = 0, then 



1. if (1 + 7) (a + 7) < 1 then almost surely X n = for n sufficiently large, and the proportion 
of isolated vertices in G n tends to 1 almost surely as n — >■ 00; 

2. if (1 + 7) (a + 7) > 1 then there is q > such that the probability that X n — >■ 00 as n — >■ 00 
is q and the probability that X n — > as n — > 00 is 1 — q. 



Proof. We note that (X n ) follows a Smith- Wilkinson branching process in random environment, 
[9]. The environmental variables which determine the random environment are the random 
variables £„, with the offspring distribution of the branching process at time n having mean 
1 + 7 if £rt+i = 1 and a + 7 if £, n +i = 0. Hence, by Theorem 3.1 of [9], the branching process 
dies out with probability 1 if \ log(l + 7) + ^ log(a + 7) < 0, i.e. if (1 + 7)(a + 7) < 1, and the 
branching process dies out with probability strictly less than 1 otherwise, hence there is positive 
probability that X n — > 00 as n — > 00. To see that the proportion of isolated vertices tends to 1 
almost surely when (1 + 7)(a + 7) < 1, note that the proportion of isolated vertices is increasing 
(as if a vertex v is isolated in G n both vO and vl are isolated in G n +i) and therefore must 
converge to some value, which cannot be less than 1 as the degree of a random vertex converges 
to zero almost surely. □ 

Proposition 13. (a) If f3 > and (1 + 7) (a + 7) < 1, the degree distribution of the graph 

(n) 

converges to the stationary distribution of the Markov chain in the sense that if we let p d 
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be the proportion of vertices in G n with degree d, and let X be a random variable with the 
stationary distribution of the Markov chain, then p^ — >■ P(X = d) as n — > oo, almost 
surely, for all d G No- 

(b) If (3 > and (1 + 7) (a + 7) > 1 then p^ — > as n — > 00, almost surely, for all d £ No- 

Proof. The graph at stage r contains 2 t vq vertices. We then consider the edges of G r+S in 
two sets: those which are between descendants of the same vertex in G r , and those which are 
between descendants of different vertices in G r . For the former, the appearance of edges between 
descendants of one given vertex is independent of what happens to the descendants of the other 
vertices, so we can model these edges of G r+S as consisting of 2 t vq independent copies of G s , 
where (G n ) ne j$ represents the process as evolved from a single vertex with no edges, Go. Then, 
by Chebyshev's inequality and a Borel-Cantelli argument, as r — > 00 proportions of vertices in 
G r+S with degree d excluding connections to descendants of different vertices in G r converge, 
almost surely, to P(X S = d). In the (1 + 7) (a + 7) > 1 case we know that P(X S = d) —> as 
s — > 00 for all d, and the actual degree of a vertex is bounded below by the degree excluding 
some connections, so this is enough to prove (b). 

To complete the proof of (a), we need to consider edges between vertices which are descendants 
of different vertices in G r . We couple the process, starting from G r , with a process with /? = 
by removing all edges between a vertex and its offspring, and all edges descended from such 
edges. The edges thus removed from G r+S will all be between vertices descended from the same 
vertex in G r , so all edges between vertices descended from different vertices in G r are present in 
the f3 = version. But by Proposition 12 the proportion of vertices in G r + S which have non-zero 
degree in the f3 = version tends to zero as s — >■ 00, and so this also applies to the proportion of 
vertices in G r + S which have edges connecting them to vertices with a different ancestor in G r . 

Hence as both r and s — > 00 the proportion which have degree d converges to P(X S = d). □ 

Finally we can put the Propositions above together to deduce Theorems 1 and 2. 

Proof of Theorem 1. Theorem 1 follows from Proposition 12; in the supercritical case where 
(1 + 7) (a + 7) > 1 the probability that a randomly chosen vertex in the graph is isolated tends 
to 1 - q < 1. 

Proof of Theorem 2. Theorem 2 follows immediately from Propositions 8, 9, 10 and 13. 
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